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ABSTRACT 



Context. The synchrotron radio morphology of bilateral supernova remnants depends on the mechanisms of particle 
acceleration and on the viewing geometry. However, unlike X-ray and 7-ray morphologies, the radio emission does not 
depend on the cut-off region of the parent electron population, making it a simpler and more straightforward tool to 
investigate the physics of cosmic ray production in SNRs. 

Aims. Our aim is to derive from the radio morphology tight constraints on the direction of the local magnetic field and 
its gradient, and on the obliquity dependence of the electron injection efficiency. 

Methods. We perform a set of 3D MHD simulations describing the expansion of a spherical SNR through a magnetized 
medium with a non-uniform ISMF. From the simulations, we derive non-thermal radio maps and compare them with 
observations of the SN1006 remnant. 

Results. We find that the radio morphology of SN1006 at 1 GHz is best-fitted by a model with quasi-parallel injection 
efficiency, a magnetic field aspect angle of 38° =b 4° with the line of sight, and a gradient of the field strength toward 
the galactic plane, higher then the expected variations of the large scale field of the Galaxy. 

Conclusions. We conclude that the radio limbs of SN1006 are polar caps, not lying in the plane of sky. The study of 
the synchrotron radio emission of SNRs is of crucial importance to derive information on the galactic magnetic field in 
the vicinity of the remnants, and more hints on the actual injection efficiency scenario. 

Key words. Acceleration of particles; Shock waves; ISM: supernova remnants 
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1. Introduction 

The radio morphology of supernova remnants (SNRs) may 
be very informative on the conditions of the magnetized 
environments in which the blast- wave expands and, more 
in particular, on the acceleration processes which occurs at 
the shock front, and which give rise to the energetic elec- 
trons ultimately responsible of the synchrotron emission in 
the radio and X-ray band. In particular, the radio emis- 
sion, unlike the X-ray emission which depends on both the 
efficiency of electron injection into the acceleration process 
(i.e. the fr action of particles injected from the thermal gas, 
as in e.g. iBlasi et aD l2005f ) and on the rapidity of accel- 
eration to high energies, is insensitive to acceleration-rate 
issues. In fact, the time required to accelerate an electron 
emitting synchrotron emission at a frequency u is about 
3 X 10-^ri(u/GHz)^-^(B/ 100uG)-^-^(Vj3. 000 km s'^)'^ 
yr (|Uchivama et al.|[2007[ ), where B is the magnetic field, Vs 
is the shock speed, rj is the ratio between the mean free path 
of the particles along the magnetic field line to the gyrora- 
dius, and the parameters are normalized to typical young 
SNR values. F or the archetypical remnant SN1006, we have 
B ^ 100 uG (IVolk et ap l2005L lBalletll2006L [Morlino et al ' 
2010",'Petruk et al.""201lllJ/5 - 4000 km s-^ fMoffett et aT 
1993. Katsuda et al. 2009|), and 7? - 1 - 10 (Petru k et al 



to SNR evolutionary timescale 1 kyr). As a result, vari- 
ations of radio morphology with obliquity point to electron 
injection physics alone, making it the best regime by which 
to investigate this issue. 

In this context, the remnants showing a radio shell 
with two opposite bright and regular limbs separated by 
two minima of emission (often name d bilateral or barrel- 
shape d or bipolar, B SNRs hereafter, Kesteven fc Caswelll 
198i iGaenslerl [19981) are considered ideal laboratories, be- 
cause their morphology is definitely not heavily affected 
by small scale inhomogeneities which may make the inter- 
pretations rather difficult. A point-like supernova explosion 
which occurs in a uniform magnetized medium with con- 
stant value and direction of the interstellar magnetic field 
(hereafter ISMF) yields a symmetric radio BSNR whose 
bright limbs are located where the magnetic field is parallel 
or perpendicular to the shock normal, if the injection ef- 
ficiency is quasi-parallel or quasi-perpendicular /isotropic^, 
respectively, and if the ISMF direction lies in the plane 



201 li) . thus, for radio emitting electrons, the resulting accel- 
eration time (less than a year) is extremely rapid compared 



^ Here, perpendicular and parallel refers to the angle between 
the shock normal and the pre-shock magnetic field (this angle is 
called the obliquity angle). In particular, the injection is called 
quasi-parallel (quasi-perpendicular) if its efficiency is maximum 
where the obliquity angle is 0° (90°, iFulbright Revnoldsl 
Il 99 0l) . The injection is called is otropic if its efficiency does not 
depend on the obliquity angle (|Ellison et al.l Il995l . IVolk et al] 
[2003 ). 
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of sky, whereas different configurations occurs at different 
aspec t angles ( Fulbright fc Revnoldsl 1 19901 : lOrlando et al.l 
I2QQ7D 

However, in the real univ erse, many SNRs evolve in 
a non-uniform ISM (see e.g. ' Hnatvk fc Petru5 1 19991 for 
an analytical treatment of SNR in nonuni form medium) 
and BSNRs are therefore often asymmetric. Orl ando et al.l 
(I2QQ7) (hereafter Paper I) have generalized the study of 
iFulbright fc Revnoldsl (|l99Qf ) to these cases by considering 
the explosions occurring in a large scale gradient either of 
density or magnetic field, showing that most of the asym- 
metries observed in true BSNRs can be recovered with this 
model. In particular, the radio morphology loses one axis 
of symmetry, and the limbs are not equally bright (if the 
gradient runs across the limbs) or they are not opposite 
and converge on the side in which either the density or the 
magnetic field is increasing (if the gradient runs between 
the limbs). 

It is clear that the morphology of BSNRs is tightly cou- 
pled with the magnetized environments in which the shock 
expands, and it is of particular interest to note here the 
preference of BSNR symmetry axifl to be oriented para l- 
lel to the galactic plane, as reported bv iGaenslerl (|l998f ). 
It seems therefore possible to study their morphology to 
derive the geometry of the surrounding magnetic field, and 
thus shedding more light on the mycrophysics of the particle 
acceleration processes which undergoes at the shock front. 
In this respect, BSNRs may be considered as a probe for 
structure of Galactic MF on scales of ^ 10 pc, much lower 
than the expected large scale variation of t he field predicted 
by global models (Prouza fc Smida 20031 iKachelrieB et al.„ 
120071 . ISun et al.ll2008F ! 

The remnant of SN1006 seems to be the object in which 
this kind of study may be more fruitful. The uniform en- 
vironment and the bright limbs visible in most of the elec- 
tromagnetic spectrum makes it a real case study in the 
field of particle acceleration me chanisms in strong shocks. 
Indeed, iRothenfiug et alj (|2004l ). using a simple and pow- 
erful geometrical argument applied to the XMM-Newton 
X-ray image of SN1006, based on the ratio between the cen- 
tral and the rim luminosity, showed that the bright limbs 
are polar caps (instead of an equatorial belt) and that, 
therefore, the magnetic field is oriented across the limbs, 
in the NE-SW direction. This argument seems to break 
the dicotomy between the two competing scenarios of the 
dependence of the injection efficiency on the obliquity an- 
gle (the angle between the magnetic field direction and 
the shock normal), preferring a situation in which the in- 
jection is most efficient when the field is along the shock 
normal (quasi-parallel scenario) over the situation in which 
the field is perp endicular (quasi-perpendicular scenario, see 
iFulbright fc"^v nolds 1990). 

On one hand, thes e findings were somehow expected 
since iVolk et al.l (|2003D already pointed out that in SN1006 
the injection should be maximum at parallel shocks. 
However, in the light of the uncertainties related to the de- 
tails of the acceleration processes, several authors still con- 
sidered the quasi-perpendicu lar scenario a viable option: 
iFulbright fc Revnoldsl (|l990f ) argued against the quasi- 
parallel scenario pointing out that quasi-parallel models 
often give rise to unobserved morphologies in the ra- 



^ The axis running parallel to the limbs in the plane of the 
sky. 



Table 1. MHD models of SN1006 used in this work 



Name 


V|B[2 


UNIFORM_B 


1.0 


GRADl 


1.2 


GRAD2 


1.4 


GRAD3 


2.0 


GRAD4 


4.7 





^ The relative variation of the magnetic field strength over 
a scale of 10 pc 

dio band (see also ' Orlando et al.l [2007f ): lYamazaki et al.l 
(2004) still considered both models to explain the ob- 
served width of sel ected fila ments of SN1006 observed by 
iBamba et aH (|2003D . whereas Amano fc Hoshinol (|2007[ ) de- 
veloped a quasi-perpendicular model which agrees very well 
with the sam e data. 

Recently, iPetruk et al.l (|2009b f) have further investi- 
gated this issue, showing that, in the framework of a sim- 
ple model of SN1006 in terms of a point-like explosion oc- 
curring in a uniform density and uniform magnetic field 
medium, there is no way to reconcile the quasi-parallel sce- 
nario and the SN1006 morphology as observed in the radio 
band. The same conc lusion was reached more recently by 
ISchneiter et al.l (|2010l ). The contradiction in the interpreta- 
tion of the radio morphology (suggesting isotropic or quasi- 
perpendicular injection efficiency scenario) and X-ray mor- 
phology (suggesting quasi-parallel scenario), and the fact 
that SN1006 shows both the kinds of asymmetries discussed 
in Paper I (converging limbs and different surface bright- 
ness), has prompted us to investigate the effects of non- 
uniformity of the environment on the obser ved properties 
of this remnant, capitalizing on the work of lOrlando et al.l 
(|2007[ ) on asymmetric BSNRs. 

In Sect. 2, we brieffy describe the MHD model used to 
reproduce the remnant of SN1006, which includes a small 
gradient of the magnetic field. In Sect. 3, we introduce the 
methodology for the comparison between modeled and ob- 
served radio images of SN1006, while in Sect. 4 we discuss 
the results of the comparison, showing that, not only the 
model with the gradient of |B| is able to reconcile, for the 
first time, the interpretation of radio and X-ray morphology 
of the remnant, but it also provides stringent constraints on 
the overall geometry of the field in the vicinity of the rem- 
nant. 



2. The model 

We adopted t he 3D MHD SNR model described by 
lOrlando et al.l (|2007[ ) to which the reader is referred for 
further details, and we use the FLASH code ( Frvxell et al.l 
.200 0) for the model implementation. Since there is accu- 
mulated evidence that the density around SN1006 is fairly 
constant, with the exception of the NW s ector w here an en- 
counter with a dense cloud is occurring (|Acero et al. 20071 
LMiceli et al. 2009), we argue that a model with a |B| gradi- 
ent is more appropriate than that with a density gradient 
to describe the asymmetries of the limbs of SN1006. Given 
the pronounced slant ness of the limbs, we conclude, fol- 
lowing Paper I, that a prominent component of V|B| is 
directed parallel to the limbs (and therefore perpendicu- 
lar to the galactic plane). The remnant was modeled as a 
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point-like explosion of 1.4 solar masses of ejecta having a 
kinetic energy of £^ = 1.3 x 10^^ erg occurri ng in a uni- 
form density medium with no = 0.05 cm~^ (|Acero et al.l 
I2007D . The ISMF has a value of 30 /iG in the environ- 
ment of the explosion site. Such a high value for the lo- 
cation of SN1006 in the Galaxy, has been chosen in order 
to take into account, as a first approximation, the effects 
of magnetic field amplification, so to have a post-shock |B| 
of the order of a few hundreds of /iG in our simulations, 
in a greement with observa tions of some SN1006 filaments 
(e.g. lYamazaki et"all 12004 ). The direction of MF is along 
the X axis while its gradient is along the Z (X) axis in case 
of quasi-parallel (quasi-perpendicular) injection efficiency 
(see. Fig. 1 of Paper I for more details). V|B| values used 
in our simulations are reported in Table [H The simulations 
were stopped at t = 1000 yr, checking that the shock veloc- 
ity and remnant radius are compatible with the observed 
values (4600 km s~^ and 9.6 pc, respectively). 

We synthesized t he radio em i ssion using the 
REMLIGHT code of lOrlando et all (|201lf ). to which 
the reader is referred for more details. In particular, the 
synchrotron radiation from relativistic electrons distributed 
with a power law spectrum N(E) (x KE~^ i s computed 
according to iGinzburg fc Svrovatskiil (|l965f ) with the 

formula ex KB^^^'^\~^ ^ where a = (5 — l)/2, and 
is the component of the magnetic field perpendicular to 
the line-of-sight (LoS). The coefficient K ex eV^ includes 
all the dependencies from th e obliquity and the shock 
velocity iVsh) as described bv lOrlando et al] (|2007f ). and, 
in particular, the injection efficiency is e oc sin^QBn2 and 
6 cx cos^QBn2 in the quasi-perpendicular and quasi-parallel 
case, where 0Bn2 is the obliquity angle between the shock 
normal and the post-shock magnetic field. This particular 
ch oice of the obliquity d e penda nce derives from the work 
of iFulbright fc Revnoldsl dlQQOD. and it is largely adopted 
in th e literature fe.g. 'Orla ndo et aDl2007l ISchneiter et al.l 
2OIOI, Orlando et al. 2011; f or alternative formulations see 



Volk et all 12003 and Petruk_elal]|2009a). The parameter 
h has been chosen equal to -1.5, reflecting a situation 
in which stronger shocks inject particle more efficiently, 
but we have verified that our results are not affected by 
ch oices of h in th e range (-1.5, +2.0), the range explored 
bv lOrlando et all (j2007t ). 



3. Comparison between models and observations 

3.1. The angles defining the viewing geometry 

Several works in the literature show that the observed mor- 
phology of a SNR emitting synchrotron emission is strongly 
affected by the angle between the direction of B and the 
line of sight, 0b- Moreover, in our circumstance, since we 
have also a gradient of the magnetic field, we are forced to 
consider the dependence of the observed morphology from 
the angle between the direction of the gradient and the 
plane of the sky, 0v|B|- The definition of 0b and 0v|b| 
are sketched in Fig. [Tl To the purpose of this work, we 
have chosen to synthesize our maps in the following way: 
the starting configuration has the MF perpendicular to the 
LoS (0B = 90°), and V|B| lying in the plane of the sky (the 
plane P in Fig.[T]) and parallel to the lobes. First, we apply 
a rotation of 0v|b| degrees around the axis in the plane of 
the sky passing through the center of the remnant and per- 
pendicular to the limbs, where positive angles means that 



LoS 




Fig. 1. The angles defining the viewing geometry, namely 
0B (the magnetic field aspect angle) and 0v|b| (the rota- 
tion angle of the gradient). P is the plane of the sky, on 
which the SN1006 are projected (sketched as light grey). 
See text for more details. 



region of increasing B are closer to us. Then, we apply a 
rotation of 90° — (j)^ degrees around the axis in the plane 
of the sky passing through the remnant center, and parallel 
to the limbs (axis A in Fig. [T]). Others rotation schemes 
gives similar results. The adopted values of 0b are from 
0° to 90° in step of 2° and the ones of 0v|b| are from 0° 
to 90° in step of 15°. Therefore, for each model, we have 
generated 315 maps encompassing all the combinations of 
the relevant angles. 



3.2. The morphological parameters 

The exploration of the parameter space is a challenging 
task, because it involves the morphological comparison be- 
tween many model maps and the true SN1006 image. We 
devised a simple procedure which involves the comparison 
of the value of 4 morphological parameters calculated both 
in the synthesized radio emission maps and on the observed 
radio map of SN1006. The parameters are: the ratio be- 
tween the maximum and the minimum of surface brightness 
{Sh) around the rim (^4); the ratio between the maximum 
around the rim and the center of the remnant (C), the ratio 
of Sh of the two bright limbs {Rmax)] the angular separation 
between the limbs maxima (Op ). The paramet e rs A, RmM,x 
and Od have been also used in lOrlando et al.l (|2007l ), and 
we refer to that paper for further discussion. The C param- 
eter is introduced here to measure the luminosity contrast 
between the brightest rim and the center. Note that the 
C para meter does not co r respo nd to the R^^/s parameter 
used by iRothenfiug et al.l (|2004l ) to exclude the equatorial 
belt scenario for SN1006, because C is measured in a small 
circular region in the radio map (see below). While Rmax 
and Od gives a measure of the asymmetries (if any) in the 
remnant morphology, A and C give the idea of the rem- 
nant being centrally filled (C > 1, A ~ 1) or barrel-shaped 
(C < 1, A > 1). This choice of evaluating the emission 
in selected regions makes the method less sensitive to the 
exact functional form of the obliquity dependance of the 
injection efficiency introduced in Sect. 2. 



4 



Bocchino et al.: Local B field around SN1006 



-41:40:00.0 



50:00.0 



-42:00:00.0 



10:00.0 




03:20.0 15:02:00.0 




0.0001 0.0002 0.0003 0.0004 0.0005 



Fig. 2. A radio map of SN1006 at 1.4 GHz, with the 
four region in which the surface brightness of the remnant 
has been averaged for the purpose of computation of the 
four morphological parameters: A = RIMmax/ Rl^mim 
C = center /RIMmax, Rmax RIMmaxi/RIMmax2 and 
Od: t he aperture a n^le between the limb maxima. Adapted 
from lPetruk et all ([2QQ9b ). 
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Fig. 3. Comparison between synthesized radio maps in the 
case of a uniform B at an age of lOOOyr (black, green 
and red dots for quasi-parallel, quasi-perpendicular and 
isotropic injection efficiency, respectively) and the observed 
radio map (black diamond, the inset show a zoom of the 
plot around this point) by means of the use of the A and C 
morphological parameters introduced in Fig. [2l Each curve 
represent the values spanned by the parameters when the 
aspect angle goes from 0° (rightmost part) to 90° (left- 
most part). 



nariofl In this case, the plot suggests that quasi-parallel 
models does not reproduce at all the observed parameters, 
unlike the quasi-perpendicular and isotropic scenarios, with 
a best- fit (/)B ^ 70°. This is i n excellent agreement with the 
work of lPetruk et al.l (j2009b[ ) , which derived (j)-B by compar- 
ing the azimuthal brightness profile with models in uniform 
magnetic field. Therefore, we conclude that our method is 
robust. 



The parameters A, C and Rmo.x are measured for 
SN1006 using the radio map of iPetruk eraP (|2009b ). by 
averaging the Sh value on circular regions of 45^'radius (Fig. 
[21). The values we got are A = 2.1 ± 0.2, C = 0.36 ± 0.03, 
Rmax = 1.2 ± 0.1, and Od = 135° ± 10°. For the measure- 
ment of the parameters in the synthesized radio maps we 
used an automated procedure to find the maximum of the 
two limbs, the minimum between limbs along the rim and 
the central position. Then we used an average in circular 
regions whose radius is the same percentage of the rem- 
nant radius used in the true radio image (5%) to minimize 
fiuctuations. Since the model images have been synthesized 
using the same number of pixels per remnant radius as the 
true radio map, this procedure ensure that the model and 
observed values of the parameters are comparable. 

3.3. Test of the procedure 

To test the procedure, we run the model with the initial pa- 
rameters discussed in Sect. 2 but assuming a uniform mag- 
netic field, in the quasi-parallel, quasi-perpendicular and 
isotropic injection efficiency scenarios. This situation is al- 
ready somehow in contradiction with the observed values 
6d = 135° and Rmax = 1-2, because V|B| = implies 
Od = 180° and Rmax = 1- However, this MF configu- 
ration has been used in the literature, and it is a good 
test case for our procedure. In Fig. [3l we show the A — C 
scatter plot for various (j)-B angles, in the uniform B case 
and quasi-parallel, quasi-perpendicular and isotropic sce- 



3.4. Results 
3.4.1. Polar caps 

In Fig. m we show the A — C, Rmax — C and Od — C scat- 
ter plots computed in the synthesized radio maps of our 
SN1006 model including a weak gradient of the magnetic 
field (GRAD2 in Table [U i.e. B varies of a factor 1.4 over 
10 pc scale), and assuming that the electron injection effi- 
ciency is maximum where the obliquity angle is 0° (quasi- 
parallel scenario). This configuration has also been referred 
to as "polar caps" . We also overplotted the values of the pa- 
rameters derived usin g the true radio map of the remnant of 
iPetruk et al.l (|2009b[ ) . Though the agreement is not exactly 
perfect in all the plots, we note that we can define a very 
limited region of the parameter space (0b,^v|b|) which is 
compatible with the observed values of A, C, Rmax and Od- 
This means that the observed radio morphology of SN1006 
is overall compatible with a quasi-parallel scenario for this 
remnant, if we include a weak gradient of B. 

A remarkable result is that the comparison between the 
quasi-parallel model and the observation strongly exclude 
a situation in which the polar caps are in the plane of 
the sky (^b = 90°) or along the line of sight (^b = 0°). 
The latter geometry would cause a centrally brightened ra- 
dio nior2holog^Y_2nstea^^ limb s, as already di s cusse d 
bv lFulbright fc Revnoldsl (|1990[ ) and lOrlando et al.l (|2007D . 
The best-fit values of the relevant angles derived from Fig. 



^ In the uniform B case, Rmax — 1 and Od — 180° by defini- 
tion. 
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Fig. 5. Top panel: Synthesized radio image at 1 GHz for 
our best fit model (GRAD2 in Table [TJ and quasi-parallel 
injection efficiency scenario, 0b = 38° and (/)v|b| = 15°), 
000 smoothed with a sigma of 2^ Bottom panel: radio map of 
SN1006 of Fig. [21 but smoothed with the same sigma used 
for the model image. 



Fig. 4. Top panel: A — C scatter plot for our model GRAD2 
(Table [T]), assuming quasi-parallel injection and an age of 
1000 yr. Different colors correspond to different values of 
0v|B| (from 0° to 90° in steps of 15°), whereas each dot 
corresponds to a given value of 0b (from 0°, right, to 90°, 
left, in steps of 2°). We overplotted the (A,C) values mea- 
sured for SN1006 (diamonds; the inset shows a zoom near 
this point). Middle panel: Same as top panel but for the 
Rmax — C pair. Lower panel: Same as top panel but for the 
9d — C pair. 



m and a conservative estimate of their uncertainties, are 
0B = 38° ± 4° and 0v|b| = 15° ± 15°. 

The synthesized radio map of the best fit model is shown 
in Fig. [5] (top panel), along with the observed radio map 
(bottom panel). Both maps have been slightly smoothed (2' 



sigma, ~ 1/15 the remnant size) to focus the comparison on 
the general obliquity trends. In fact, we are not interested 
in reproducing the fine structure of the thermal emission 
from the ejecta which is visible in the high-resolution radio 
observation. We note that the large scale structures of the 
observed radio emission are very well recovered by the best- 
fit model and the two images looks similar. 

The comparison with other gradients, namely the 
GRADl (weaker gradient), GRAD3 and GRAD4 models 
(stronger gradient), shows that it is not possible to find 
a satisfactory fit for all the parameters. In particular, for 
higher values of the gradient, the angle is always smaller 
than observed in most of the cases, and for lower values is 
always larger, so GRAD2 is the only model which gives us 
an overall good fit. The variation of the MF in SN1006 can 
therefore be bracketed between 1.2 and 2.0 (see Table [1]). 
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3.4.2. Equatorial belt 

In this section, we explore the same MHD model as in the 
previous section, but now assuming that the electron in- 
jection efficiency is maximum where the obliquity angle is 
90° (quasi-perpendicular scenario). This configuration has 
also been refereed to as "equatorial belt" . We produced the 
synthetic images and we repeated the analysis described in 
the previous paragraphs. In this case, the direction of the 
V|B| is aligned with the direction of £0, so the angle ^v|b| 
is the same as (/)b, and we did not considered it any fur- 
ther. Another consequence of this is that Rmax is always 
1 (in disagreement with observations), so this parameter 
cannot give any diagnostics. The results are shown in Fig. 
ini and it seems that a good fit can be found for (j)-B ^ 70° 
(a value in agreeme nt with the azimuthal profile analysis of 
iPetruk et al.ll2009b[ ). even if the model points in the Od — C 
diagram seems to be more distant to the observed points 
than any scatter plot in Fig. [H thus indicating a better fit 
in the quasi-parallel case then in the quasi-perpendicular 
one. 



4. Discussion 

In this work, we introduced a new method of comparison 
between models and observations of radio SNRs shells that 
can provide useful information on the 3D structure of the 
remnant and on the unperturbed ambient magnetic field in 
which the remnant expands. The method is based on the 
calculation of 4 morphological parameters. In the case of 
SN1006, we found that the method gives a good fit if the 
injection efficiency is quasi parallel and the ambient mag- 
netic field is characterized by a gradient of its strength. In 
addition, the method allowed us to constrain the viewing 
geometry angles of th e magnetic field (0b_= 38° ± 4° and 
</>v|B| = 15° ± 15°). ISchneiter erall (|2010[ ) showed that 
the quasi-perpendicular model is a better fit to the R^^/s 
value (the ratio betwee n the inn er region and the limb re- 
gion flux introduced by iRothenfiu g et al. 2004) of SN1006 
in the radio band, but they did not properly considered as- 
pect angles ^ 90° and the non-uniform magnetic field. We 
verified that our model gives R^^/s = 0.16 for 0b = 90°, 
but R^/s = 1.00 for <p-B = 38°, in agreement with the value 
measured in the new radio map of Fig. [2] (i^^yg = 1.01). 

The symmetry axis of SN1006 is perpendicular to the 
galactic plan^E, its distance is 2.2 kpc ([Winkler et al.|[2QQ3[ ) 
and its galactic latitude is 14.6°. Combining all this pieces 
of information, we can plot the direction of the magnetic 
field and its gradient in a 3D representation of the Galactic 
disk. This is shown in Fig. [71 

Quite remarkably, our best-fit values of 0v|b| and 0b 
imply that the direction of V|B| resulting from the fit to 
the quasi-parallel model points down toward the Galactic 
plane and the direction of B is aligned with the direction of 
the spiral arm near the remnant. This is indeed a very rea- 

^ This is strictly necessary because ' Orlando et al.l (|2007l l has 
shown that the gradient increasing direction is parallel to the 
limb, so quasi-perpendicular injection efficiency implies B per- 
pendicular to shock normal at the limbs and therefore B parallel 
to V|B|, 

^ This makes SN1006 rather peculiar, since the majority of 
bilater al SNRs have t heir symmetry axis parallel to the galactic 
plane (|Gaensledll9"98h . 




0,01 OJO 

C (center/rim^^J 

Fig. 6. Top panel: Same as Fig. H] (top panel) but for the 
GRAD2 model and quasi-perpendicular scenario. There is 
no ^v|B| angle shown in this case, because B is always 
aligned with V|B|. Lower panel: Scatter plot oi Od — C 
parameters. We overplotted the values observed in SN1006 
(diamonds). 




Fig. 7. Artist view of the magnetic field at the location 
of SN1006 in our Galaxy. The red and green arrows mark 
the directions of the field and its gradient as derived by 
the best-fit model for the synchrotron radio emission of 
SN1006, assuming a quasi-parallel scenario, derived in this 
work. The figure is in scale. 



sonable configuration for the magnetic field of the Galaxy, 
so it is tempting to conclude that we are sampling the large 
scale field of the Galaxy. However, we note that the best-fit 
value of V|B| (i.e. a variation of a factor 1.4 over 10 pc, 
model GRAD2 in Table [1]) seems to be high if compared 
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with the values predicted by current models of galactic MF. 
We have used the model of'S un et al.l (|2QQ8[ ) to compute the 
expected variation of the azimuthal field in the Galaxy at 
the location of SN1006. In particular, we have used their 
ASS+RING model which seems to be one of the favored by 
rotation measures of pulsars in the plane. At the SN1006 
distance of ~ 550 pc from the plane, the large-scale field is 
expected to vary of a factor of 1.4 on scales of more then 
100 pc, much larger than the scale- length required by our 
lowest gradient model (GRADl in Table [1]). This opens 
up the possibility that we are a ctually sampling the ran - 
dom magnetic field component. Minter & Spangler (1996) 
reports length scales for this component of the order of few 
pc, which is compatible with the variations we derive. 

5. Summary and conclusion 

The synchrotron radio emission of the archetypical bilat- 
eral supernova remnant SN1006 is compared against a 
MHD model for this remnant, including a gradient of the 
magnetic field strength, particle acceleration and its obliq- 
uity dependent characteristic, namely quasi-perpendicular, 
quasi-parallel and isotropic scenarios for the injection ef- 
ficiency and a rough treatment of magnetic field amplifi- 
cation. In order to explore the parameters space, which is 
very large due to the necessity to include the viewing ge- 
ometry in the model-data comparison, we developed a sim- 
plified procedure based on the computation of four mor- 
phological parameters. We have found a very good fit with 
a model assuming quasi-parallel injection efficiency and a 
configuration of the ambient magnetic field characterized 
by a variation of its strength of abut a factor of 1.4 (1.2- 
2.0) over a scale of 10 pc, and by viewing geometry angles 
(/)B = 38° ± 4° (between the direction of B and the line 
of sight), and ^v|B| = 15° ± 15° (between the plane of 
the sky and the direction of the gradient of the magnetic 
field). A worse fit is obtained in the quasi-perpendicular sce- 
nario. The overall morphology of the observed radio emis- 
sion at 1.4 GHz is correctly recovered by our best-fit model. 
Therefore, we conclude that the SN1006 limbs are polar 
caps, significantly tilted with respect to the plane of the 
sky. The projected direction of B and V|B| in the Galaxy 
are along the spiral arm and toward the plane respectively, 
which is in very good agreement with the expected direction 
of the large scale galactic B. However, the implied gradient 
is too high to be associated to the large scale galactic B 
at the location of the SN1006 remnant and more typical of 
the random magnetic field components. 

The application of our method to selected samples of 
bilateral supernova remnants may yield independent esti- 
mates of the geometry of the galactic field at several loca- 
tions, which can be useful to understand the field topology 
in our Galaxy. 
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